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Physical characteristics of composite materials have been 
investigated with increasing frequency in response to the rise 
in industry requirements for materials with specific properties. 
They include work on thermal and electrical conductivity of 
polymeric materials containing powders of different conduct- 
ing fillers as the dispersed phase (Cullen et al., 1988; Holrook, 
1988; Dinwiddie and Onn, 1989; Rule et al., 1990). 

The conductivity of such composite materials has essentially 
a nonadditive characteristic depending to a great extent on the 
distribution of the filler particles in the matrix, that is, de- 
pending on the composition structure. In a real composite, 
single, isolated particles, complexes of associated particles, and 
chain structures may exist in some combinations (Gul and 
Shenphil, 1984). 

The aim of this work is to show the construction of a model 
that offers the possibility to describe such systems, considering 
the particle fractions existing in three structural phases. Several 
theoretical models were suggested recently that determined the 
relation between microstructure and effective conductivity of 
filled composites (Beasley and Torquato, 1986; Smith and Tor- 
quato, 1989; Lu and Kim, 1990). However, simplifLed semi- 
theoretical models, especially the Lewis and Nielson (Nielson, 
1973) method including structural parameters, are used often 
by experimentalists (Hoffman et al., 1989) for predictions of 
thermal conductivity of composite materials. To describe their 
experimental data, Kusu and Cornelinssen (1975) suggested 
using the formula for the system of plates parallel to the flow 
with the introduction of an empirical coefficient. There seems 
to be no basis for using such a simplified model, which reflects 
inadequately the real, three-dimensional composite structure. 

According to the studies of de Vriese (1955) and Rivkin 
(1989), the Bruggemann model (1935) proved to be better for 
the calculation of the thermal conductivity of filled polymers 
for isolated, uniformly-distributed inclusions, while that of 
Dul'nev (see DuFnev and Novikov, 1977) was better for mu- 
tually-penetrating components. The models, however, describe 
only the conductivity of limiting structures without offering . 
the possibility of going into intermediate variants. 
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Description of the IModel 

Take the heterogeneous system, in which the particles may 
exist in three different phases: as isolated inclusions, in isolated 
aggregates of associated particles, and in chain structures (skel- 
eton). 

The relationship of phases is characterized by parameters a, 
and Ug, which are defined by particle fractions in the skeleton 
and the aggregates, respectively. Thus, the volumetric fractions 
of the skeleton and the aggregates Va are the tightly packed 
systems of the filler particles expressed as follows: 



The heterogeneous system is considered a structure with 
mutually penetrating components: one the conducting skeleton 
and the other the medium, containing isolated particles and 
the aggregates. The thermal conductivity calculation of such 
a system is carried out in four stages. 

1 . The thermal conductivity calculation of skeleton and 
aggregates \, 

2. The effective thermal conductivity calculation of the 
polymer base with isolated inclusions (aggregates and skeleton 
are not taken into account). 

3. The effective thermal conductivity calculation of the 
system of isolated aggregates in the medium with thermal con- 
ductivity \m. 

4. The effective thermal conductivity calculation of the sys- 
tem with mutually-penetrating components X, which is made 
up of the conducting skeleton and the medium with thermal 
conductivity X^^. 

Figure 1 shows the order in which the structure's levels were 
considered when calculating the thermal conductivity of the 
system. The effective thermal conductivity of the aggregates 
\ and the skeleton X^, as well as that of systems containing 
the isolated inclusions X^^ and the isolated aggregates are 
given by the Bruggemann (1935) model: 
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Figure 1. Order of levels of the structure. 

The choice of formula 1, which was recommended by 
de Vriese (1955) particularly for the thermal conductivity cal- 
culation of polymer systems containing isolated, evenly-dis- 
tributed inclusions to high concentration, is based on the 
following. 

1. In the extensive experimental study of Baranovsky (1970) 
where measures were incorporated to prevent the association 
of inclusions in the polymeric matrix, the data obtained for 
the composite thermal conductivity correlate well with the cal- 
culation according to formula 1 for different fillers. 

2. Formula! was successfully applied by Pugatch and Rivkin 
(1988) for filled polymeric materials and for the calculation 
of thermal conductivity of the skeleton consisting of associated 
particles. 

When calculating the aggregate thermal conductivity in 
the formula 1 , the values were substituted: 

Xi = X,; X„=X2; i^u = t^. 
On calculating the thermal conductivity of the skeleton X,: 

Xi = Xi; Xn = X2; Vu = V2' 

The values of the particle concentration in the skeleton vi and 
in the aggregates depend on the density of their packing. 
In accordance with the recommendation by Rivkin (1989), the 
following values, v 2=1^- 0,6, were used in calculation. 
The effective thermal conductivity \m is calculated with: 

V2 



The effective thermal conductivity %m is calculated with: 



Xi = Xji^; X|i — X^; Vu — 



At the fourth stage of calculation, the Durnev model is used 
for the system with mutually-penetrating components. The 
thermal conductivity of such a system is found by the formula: 



X = \,[c^ + ^(l-c)2 + 2A:c(l-c)(A:c+l-c)-'] 



(2) 
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Figure 2. Calculated reduced thermal conductivity X/Xi 
vs. filler concentration for composites. 

They also include isolated inclusions, as wcU as (a) particle ag- 
gregates and (b) skeleton. 

k = WX, ^ ^ 

c = positive real root of 2r - 3<r + 1 = u 

which can be determined by the formula: 

c= 0.5 - sin[l/3arc sin(2M - 1)] 

where 

Figure 2 shows calculated dependencies of the reduced ther- 
mal conductivity of the composition X/X| on the filler con- 
centration for the systems containing isolated aggregates (Figure 
2a) and conducting skeleton (Figure 2b). 

A rise in the reduced thermal conductivity with an increase 
in the structural parameters of results from the growing un- 
evenness of the filler distribution. Indicated parameters in our 
model are in accordance with the two levels of the structure 
regularity: corresponds to isolated aggregates (particle ac- 
cumulations evenly distributed in the matrix); and a, corre- 
sponds to the chain structures (particle accumulations unevenly 
distributed in the matrix as the skeleton). 

In the two limiting cases of using such a model, the calcu- 
lation leads to Bruggemann's formula (0^ = 0^= 0) with no filler 
structurization (fully-uniform system) and to the calculation 
according to Dul'nev's formula (afl=0, 0^= 1), where all the 
particles pass on to the conducting skeleton (a maximum un- 
even structure of the system). 

When the average volumetric concentration of the filler in 
the composite approaches the particle concentration in the 
aggregates and the skeleton, the curves in Figures 2a and 2b 
corresponding to different a« and a, come nearer; and when 
t;js= they are combined to one point. According to the 

model, such a composite is the tightly-packed system of the 
particles evenly-distributed in the matrix. 

Discussion 

In real composites, the content of the filler in any structural 
phase depends on many factors, especially on the technology 
of the composite preparation. The process of structurization 
and the degree of its completeness are affected by the viscosity 
of the mixture, the procedure of dispersion, and the regime 
of composite hardening. There are ways to regulate the dis- 



1740 



November 1991 Vol. 37, No. 11 



AIChE JouniHl 



^1 




Figure 3. Reduced thermal conductivity of composites 
X/Xt vs. filler concentration 1/21 X2/Xts500, 

Calculation with following values of structural parameters; (a) 
0; (b) Ua- 1, af=0; and (c) a,=0. a,= 1. 

tribution characteristics of the filier particles, such as the press- 
ing of the polymer and filler powders of different dispersivity , 
which results in the formation of the conducting skeleton at 
the borders of the polymer grains (Gul and Shenphil, 1984), 
the filling of powder coatings (Pugatch et al., 1988), and the 
introduction in the polymeric matrix of the substance pre- 
venting the particle association (Baranovsky, 1970). 

In the case of composites obtained by the dispersion of the 
filler particles in a polymer melt, the character of the struc- 
turization depends on the relationship among interaction ener- 
gies ei„ e,2, and 622- Gul and Shenphil (1984) suggested an 
approach to predicting the structure that depends on the in- 
teraction energies. If the conditions en>ei2 and ^22>^i2 are 
fulfilled, a large number of aggregates are formed in composite 
which prevents formation of a continuous structure. If > eu 
and ei2>e22, the particles wiU be separated by considerable 
polymer layers and will not be inclined to associate. Conducting 
chain structures can be formed when on a part of the particle 
surface the conditions of fij2 > en are observed at the same time 
that there are sections on the surface for which the condition 
e22>ei2 is true. 

With real composites where there are, as a rule, different 
relations to the interaction energies, the structural modifica- 
tions mentioned previously can exist in some combinations. 
In most cases, when it is impossible to determine the structure 
of the filler distribution in the composite, that is, when and 
a, are unknown, we can point out the region of possible values 
of thermal conductivity. Such a region is shaded in Figure 3 
for the relation of the filler and the matrix thermal conduc- 
tivities X2/X1 = 500, Curve a in this diagram corresponds to the 
uniform distribution of the filler a«=a,=0. Curve c corre- 
sponds to the limiting, uneven distribution of the filler aa = 0 
and oi^^l. Curve b corresponds to the location of the filler in 
isolated aggregates afl= 1 and a,= 0, and takes an intermediate 
position. 

For a description of concrete experimental data, we take 
a,=a and =0. In this case, the location of the composition 
in the indicated region is determined by the volumetric con- 
centration of the filler V2 and the parameter a. In this case, 
the parameter a characterizes irregularity of the filler location 




Figure 4. X vs. for epoxy composites [6]. 

Curves a-c were calculated according to our model with the fol- 
lowing values of a,: (a) 0.2, (b) 0.4. (c) 0.6. (d) 0.8, and (c) I, 
while ota^O (symbols as in Figure 5). 



or the degree of structurization without any difference in which 
the structural phase the filler is present. To divide clearly the 
compositions into degrees of structurization, it is convenient 
to pass on to dimensionless thermal conductivity: 

(3) 



where 

X' = thermal conductivity of the umform system calculated accord- 
ing to our model with the parameters ag^a^^O 

X* = thermal conductivity of the limiting, uneven system where the 
whole of the fiUer is located in the conducting skeleton (of«=0. 

In the calculation we assume that the parameter a does not 
depend on the concentration of the filler and is constant for 
the given composition. This means that an increase in filler 
content in the polymeric matrix does not affect the character 
of the structurization. When a is constant there are no per- 
colation effects, which is confirmed by experimental data for 
filled polymers and can be explained by the presence of a 
polymer film between the p£irticles. 

Figure 4 shows the grid of X, Vj calculated behavior as well 
as experimental data obtained by Bivkin (1989) for epoxy com- 
posites. Such data representation made it possible to classify 
studied fillers in terms of structural parameter a. The results 
are shown in Table 1. Presented values of a can be used in 
the first approximation for the calculation of thermal con- 
ductivity of the polymeric composites obtained by means of 
dispersion in a melt. The experimental thermal conductivity 
data (Rivkin, 1989) and the results calculated with our model 
are shown in Figure 5 in terms of X/X, and V2. According to 
Fig. 5 and Table 1 , it is quite obvious that deviation of the 
experimental data from the model calculations are within 5- 
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Figure 5. Reduced thermal conductivity of epoxy com- 
posites X/Xi as a function of filler concentra- 
tion (7'=' 350 K). 



(1) A1 

(2) Zn 

(3) Ti 

(4) W 

(5) Ni 

{€) bronze flakes 

(7) aluminum flakes 

(8) colloidal graphite 

(9) graphite (30-50 litn) 



(tl) AlA^l (20-30) 

(12) AlA #2 (0.5-1) 

(13) Si^N, 

(14) AJN#1 (I-IO) 

(15) AIN #2 (0.01-0.1) 

(16) AIN (70% #1 + 30% #2) 

(17) B4C 

(18) SiC 

(19) ZnO 



(10) CrjOj 

Curves were calculated according to our model with \, = 0,22 
(W/m K) and the following: 

(a) X2 = 36 WAm K). a„ = 0, a,= 0.2 

(b) Xj^lOO W/(m K). a^=0, a,«0.4 

(c) X2=17 W/(m K), aa^0,as=\.O 

(d) X3=150 W/(m K), a. = 0, a,= 1.0 



10%, that is, it does not exceed experimental errors. For brev- 
ity, analytical predictions of our model will not be compared 
with extensive experimental data acquired in other works. 

As soon as parameter a characterizes the structure of the 
composite, it can be used to predict some properties connected 
with the structure. It is also possible in some cases to solve the 
reverse problem. 



Table 1, Faiers Classified in Terms of Structural Parameter 





Value of 


Maximum 




Deviation 




Structural 








Parameter 






Fillers 


a 


K 


MAX 


ZnO, C (Colloidal Graphite), 




0.10 




W, Ni, Ti, Aluminum 


1 




and Bronze Flakes 








AI2O3, CfiOj, Al. Zn. 


0.2 


0.05 




B4C, SiC 






Si3N4, AIN of Different 




0.05 




Dispersivities, 


0.4 




C (30-50 fim) 
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Notation 



c. k, u 

en 

en 



dimensionless parameters 

interaction energy of polymer macromolecules 

interaction energy of polymer and filler particles 

interaction energy of filler particles 

volumetric concentration of considered inclusions (paru- 

cles or aggregates) 

volumetric concentration of filler particles 
panicle concentrations in aggregate and skeleton 
volumetric fractions of aggregates and skeleton in the com- 
posite 



Greek letters 

structural parameter 

particle fractions in aggregates and skeleton 
effective thermal conductivity of composite 
thermal conductivity of considered matrix 
thermal conductivity of considered inclusions' 
thermal conductivities of continuous phase (polymer) and 
filler particles 

effective thermal conductivities of aggregates and skeleton 
calculated and experimental thermal conductivities of com- 
posites , 
effective thermal conductivity of the polymer base with 

isolated inclusions 

effective thermal conductivity of the system of isolated 
aggregates in the medium with thermal conductivity \m 
X ' = calculated value of X with = a, = 0 
X' = calculated value of X with a,=0, a,= I 
X = dimensionless thermal conductivity defined by Eq. 3 



a 
X 

X„ 
X„X, 
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